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Abstract. In this paper, we study the underlying geometry in the classical Hamilton-Jacobi equation. The proposed 
formalism is also valid for nonholonomic systems. We first introduce the essential geometric ingredients: a vector 
bundle, a linear almost Poisson structure and a Hamiltonian function, both on the dual bundle (a Hamiltonian system). 
Prom them, it is possible to formulate the Hamilton-Jacobi equation, obtaining as a particular case, the classical theory. 
The main application in this paper is to nonholonomic mechanical systems. Por it, we first construct the linear almost 
Poisson structure on the dual space of the vector bundle of admissible directions, and then, apply the Hamilton-Jacobi 
theorem. Another important fact in our paper is the use of the orbit theorem to symplify the Hamilton-Jacobi equation, 
the introduction of the notion of morphisms preserving the Hamiltonian system; indeed, this concept will be very useful 
to treat with reduction procedures for systems with symmetries. Several detailed examples are given to illustrate the 
utility of these new developments. 
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1. Introduction 

The standard Hamilton- Jacobi equation is the first-order, non-hnear partial differential equation, 

for a function S{t, q^) (called the principal function) and where H is the Hamiltonian function of the 
system. Taking S{t, q^) = W{q^) — tE, where E is a constant, we rewrite the previous equations as 

. dW 

^(^^^) = ^' (1-2) 



where W is called the characteristic function. Equations (1.1) and (1.2) are indistinctly referred as 
the Hamilton- Jacobi equation (see [Ullll; see also [Gj for a recent geometrical approach). 

The motivation of the present paper is to extend this theory for the case of nonholonomic me- 
chanical systems, that is, those mechanical systems subject to linear constraints on the velocities. In 



Remark |5.11| of our paper, we carefully summarize previous approaches to this subject. These tried 
to adapt the standard Hamilton- Jacobi equations for systems without constraints to the nonholo- 
nomic setting. But for nonholonomic mechanics is necessary to take into account that the dynamics 
is obtained from an almost Poisson bracket, that is, a bracket not satisfying the Jacobi identity. In 
this direction, in a recent paper [20j, the authors have developed a new approach which permits to 
extend the Hamilton- Jacobi equation to nonholonomic mechanical systems. However, the expression 
of the corresponding Hamilton- Jacobi equation is far from the standard Hamilton-Jacobi equation 
for unconstrained systems. This fact has motivated the present discussion since it was necessary 
to understand the underlying geometric structure in the proposed Hamilton-Jacobi equation for 
nonholonomic systems. 

To go further in this direction, we need a new framework, which captures the non-Hamiltonian 
essence of a nonholonomic problem. Thus, we have considered a more general minimal "Hamiltonian" 
scenario. The starting point is a vector bundle td '■ D — > Q such that its dual vector bundle 
td* '■ D* — > Q is equipped with a linear almost Poisson bracket {■, ■}£)*, that is, a linear bracket 
satisfying all the usual properties of a Poisson bracket except the Jacobi identity. The existence of 
such bracket is equivalent to the existence of an skew-symmetric algebroid structure (|-, ■]d,Pd) on 
Td '■ D — > Q (i.e. a Lie algebroid structure eliminating the integrability property), or even, the 
existence of an almost differential on '■ D — > Q, that is, an operator which acts on the 
"forms" on D and it satisfies all the properties of an standard differential except that (d^)^ is not, in 
general, zero. We remark that skew-symmetric algebroid structures are almost Lie structures in the 
terminology of jSl] (see also |35]) and that the one-to-one correspondence between skew-symmetric 
algebroids and almost differentials was obtained in [34j- We also note that an skew-symmetric 
algebroid also is called a pre-Lie algebroid in the terminology introduced in some papers (see, for 
instance, flbl \T6\ 123] ) and the relation between linear almost Poisson brackets and skew-symmetric 
algebroid structures was discussed in these papers (see also [121 113] for some applications to Classical 
Mechanics) . 

In this framework, a Hamiltonian system is given by a Hamiltonian function h : D* — > M. The 
dynamics is provided by the corresponding Hamiltonian vector field J-C^^* ( (/) = {f,h}D*, 
for all real function / on D*). Here, A^* is the almost Poisson tensor field defined from {■,-}d*- 
The reader can immediately recognize that we are extending the standard model, where D = TQ, 
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D* = T*Q and [■, -Jd is the usual Lie bracket of vector fields which is related with the canonical 
Poisson bracket on T*Q, so that is just, in this case, the usual exterior differential. Another 
important fact is the introduction of the notion of morphisms preserving the Hamiltonian system; 
indeed, this concept will be very useful to treat with reduction procedures for nonholonomic systems 
with symmetries. We remark that this type of procedures were intensively discussed in the seminal 
paper [B] by Bloch et al. 

In the above framework we can prove the main result of our paper: Theorem In this theorem, 
we obtain the Hamilton-Jacobi equation whose expression seems a natural extension of the classical 
Hamilton- Jacobi equation for unconstrained systems, as appears, for instance, in pQ. Moreover, our 
construction is preserved under the natural morphims of the theory. This fact is proved in Theorem 

Furthermore, using the orbit theorem (see [2]), we will show that the classical form of the Hamilton- 
Jacobi equation: H o a = constant, with a : Q ^ D* satisfying d^a = 0, remains valid for a 
special class of nonholonomic mechanical systems: those satisfying the condition of being completely 
nonholonomic. See Section [s] for more details and also the paper by Ohsawa and Bloch [Slj for the 
particular case when D is a distribution on Q. 

The above theorems are applied to the theory of mechanical systems subjected to linear nonholo- 
nomic constraints on a Lie algebroid A. The ingredients of this theory are a Lie algebroid ta '■ A —>■ Q 
over a manifold Q, a Lagrangian function L : A ^ M. of mechanical type, and a vector subbundle 
td : D ^ Q of A. The total space D of this vector subbundle is the constraint submanifold (see 
[7]). Then, using the corresponding linear Poisson structure on A*, one may introduce a linear al- 
most Poisson bracket on D*, the so-called nonholonomic bracket. A linear almost Poisson bracket 
on D which is isomorphic to the nonholonomic bracket was considered in [7j; however, it should 
be remarked that our formalism simplifies very much the procedure to obtain it. Using all these 
ingredients one can apply the general procedure (Theorems 4.1| and 4.12) to obtain new and inter- 
esting results. We also remark that the main part of the relevant information for developing the 
Hamilton-Jacobi equation for the nonholonomic system (L, D) is contained in the vector subbundle 
D or, equivalently, in its dual D* (see Theorems 4.1 and 4.12). Then, the computational cost is lower 
than in previous approximations to the theory. 

In the particular case when A is the standard Lie algebroid ttq : TQ —>■ Q then the constraint 
subbundle is a distribution D on Q. The linear almost Poisson bracket on D* is provided by the 
classical nonholonomic bracket (which is usually induced from the canonical Poisson bracket on T*Q), 
clarifying previous constructions dS], EH]- In addition, as a consequence, we recover some of the 
results obtained in [20] about the Hamilton-Jacobi equation for nonholonomic mechanical systems 



(see Corollary 5.9). Moreover, we apply our results to an explicit example: the two-wheeled carriage. 
On the other hand, if our Lagrangian system on an arbitrary Lie algebroid A is unconstrained 
(that is, the constraint subbundle D = A) then, using our general theory, we recover some results 
on the Hamilton-Jacobi equation for Lie algebroids (see Corollary 5.1) which were proved in [25] . 
Furthermore, if A is the standard Lie algebroid t^q : TQ Q then we directly deduce some 
well-known facts about the classical Hamilton-Jacobi equation (see Corollary 5.2). 

Another interesting application is discussed; the particular case when the Lie algebroid is the 
Atiyah algebroid : A = TQ/G ^ Q = Q/G associated with a principal G-bundle F : Q ^ 
Q = Q/G. In such a case, we have a Lagrangian function Z : A — > M of mechanical type and a 
constraint subbundle Tjj : D Q of : A = TQ/G Q. This nonholonomic system is precisely 
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the reduction, in the sense of Theorem 4.12, of a constrained system {L,D) on the standard Lie 



algebroid ta ■ A = TQ — > Q. In fact, using Theorem |4.12[ we deduce that the solutions of the 
Hamilton- Jacobi equations for both systems are related in a natural way by projection. We also 
characterize the nonholonomic bracket on D*. All these results are applied to a very interesting 
example: the snakeboard. In this example, an explicit expression of the reduced nonholonomic 
bracket is found; moreover, the Hamilton- Jacobi equations are proposed and it is shown the utility 
of our framework to integrate the equations of motion. 

We expect that the results of this paper will be useful for analytical integration of many difficult 
systems (see, as an example, the detailed study of the snakeboard in this paper and the examples 
in [21]) and the key for the construction of geometric integrators based on the Hamilton- Jacobi 
equation (see, for instance. Chapter VI in and references therein for the particular case of 
standard nonholonomic mechanical systems). 

The structure of the paper is as follows. In Section [2| the relation between linear almost Poisson 
structures on a vector bundle, skew-symmetric algebroids and almost differentials is obtained. In 
Section [3} we introduce the notion of a completely nonholonomic skew-symmetric algebroid and 
we prove that on an algebroid D of this kind with connected base Q the space H^{d^) = {/ G 
C°^{Q)/d^f = 0} is isomorphic to M. We also prove that on an arbitrary skew-symmetric algebroid 
D the condition d^f = implies that / is constant on the leaves of a certain generalized foliation 



see Theorem 3.4). For this purpose, we will use the orbit theorem. In Section |4| we consider 



Hamiltonian systems associated with a linear almost Poisson structure on the dual bundle D* to a 
vector bundle and a Hamiltonian function on D*. Then, the Hamilton- Jacobi equation is proposed 
in this setting. Moreover, using the results of Section [3} we obtain an interesting expression of this 
equation. In Section [5} we apply the previous results to nonholonomic mechanical systems and, in 
particular, to some explicit examples. Moreover, we review in this section some previous approaches 
to the topic. We conclude our paper with the future lines of work and an appendix with the proof 
of some technical results. 

2. Linear almost Poisson structures, skew-symmetric algebroids and almost 

differentials 

Most of the results contained in this section are well-known in the literature (see [TH [T5| [T6 t [Ml [35] ) . 
However, to make the paper more self-contained, we will include their proofs. 

Let td '■ D ^ Q he a vector bundle of rank n over a manifold Q of dimension m. Denote by D* 
the dual vector bundle to D and by td* : D* ^ Q the corresponding vector bundle projection. 

Definition 2.1. A linear almost Poisson structure on D* is a bracket of functions 

{-, -ji). : C°°(D*) X C°°(D*) ^ C°°(D*) 

such that: 

(i) {■,-}d* is skew- symmetric, that is, 

{^,i;}D* = -{i;,^}D^, for^,^eC^{D*). 

(ii) {■,-}d* satisfies the Leibniz rule, that is, 

{ipip',^}D' = ip{ip',ip}D' +^'W,i^}D', forip,ip',^ e C^iD*). 
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(iii) {■, ■}d* is linear, that is, if (f and ip are linear functions on D* then {(f,ip}D* is also a linear 
function. 

If, in addition, the bracket satisfies the Jacobi identity then we have that {■,-}d* is a linear 
Poisson structure on D* . 



Properties (i) and (ii) in Definition 2.1 imply that there exists a 2- vector Ad* on D* such that 

An^idip^d^P) = {ip,^P}D*, for ip,^P e C^iD*). 

Ajj* is called the linear almost Poisson 2-vector associated with the linear almost Poisson 
structure {-, -jo.. 

Note that there exists a one-to-one correspondence between the space T^td) of sections of the 
vector bundle '■ D ^ Q and the space of linear functions on D*. In fact, if X e ^{td) then the 
corresponding linear function X on D* is given by 

X{a) = a{X{TD*{a))), for a G D*. 

Proposition 2.2. Let {-, - jo. be a linear almost Poisson structure on D* . 

(i) If X is a section of td : D Q and f is a real C°° -function on Q then the bracket {X, f o 
td* }d* is a basic function with respect to the projection td* . 

(ii) /// and g are real -functions on Q then 

{/ o ■rD',go td*}d' = 0. 

Proof. Let Y be an arbitrary section oi : D ^ Q . 
Using Definition 2A_, we have that 

{X, (/ o rD*)Y}D* = (/ o TD^){X, Y}d^ + ({X, / o rD^}D*)Y 

is a linear function on D* . Thus, since (/ o to*){X,Y}£,* also is a linear function, it follows that 
{X, / o T/).}/). is a basic function with respect to td*. This proves (i). 
On the other hand, using (i) and Definition 2.1[ we deduce that 

{(/ o 'TD')Y,go Tr)*}D' = (/ o rr)'){Y,go td*}d* + ({/ o ro^^go Tr)*}D')Y 

is a basic function with respect to td*. Therefore, as (/ o td*){Y ,g o T£,*}d* also is a basic function 
with respect to r^j*, we conclude that {/ o TD*,g o te,*}d* = 0. This proves (ii). □ 

If (g*) are local coordinates on an open subset U of Q and {Xq,} is a basis of sections of the vector 
bundle t^^{U) U then we have the corresponding local coordinates {q^,Pa) on D*. Moreover, from 
Proposition |2.2[ it follows that 

{Pa,Pl3}D* = -ClpP^, {q\Pa]D-' = P^, {q\q^}D* = 0, 

with Clfs and pi real C°°-functions on U. 

Consequently, the linear almost Poisson 2-vector A^* has the following local expression 
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Definition 2.3. An skew- symmetric algebroid structure on the vector bundle td : D ^ Q is 
a M-linear bracket -Jd '■ ^(td) x ^{td) — * ^(td) on the space T{td) and a vector bundle morphism 
Pd : D ^ TQ, the anchor map, such that: 

(i) [■, ■\d is skew-symmetric, that is, 

lX,Yln = -lY,XjD, forX,Y eTiTo). 

(ii) If we also denote by pa : T^td) ^{Q) the morphism of C°°{Q) -modules induced by the 
anchor map then 

lX,fYlD = flX,YjD + PD{X){f)Y, forX,YeT{D) andfeC'^iQ). 

If the bracket |-, satisfies the Jacobi identity, we have that the pair ([-, -jDyPD) is a Lie alge- 
broid structure on the vector bundle To '■ D ^ Q . 

Remark 2.4. If {D, |-, -J^i, po) is a Lie algebroid over Q, we may consider the generalized distribution 
D whose characteristic space at a point g G Q is given by D{q) = poiDq), where Dg is the fibre 
of D over q. The distribution D is finitely generated and involutive. Thus, D defines a generalized 
foliation on Q in the sense of Sussmann D is the Lie algebroid foliation on Q associated 
with D. o 

Now, we will denote by LAy^D*) (respectively, the set of linear almost Poisson structures 

(respectively, linear Poisson structures) on D*. Denote also by SSA{D) (respectively, LA{D)) the 
set of skew-symmetric algebroid (respectively. Lie algebroid) structures on the vector bundle : 
D —>■ Q. Then, we will see in the next theorem that there exists a one-to-one correspondence 
between HAy^D*) (respectively, £?(/}*)) and the set of skew-symmetric algebroid (respectively. Lie 
algebroid) structures on td : D ^ Q. 

Theorem 2.5. There exists a one-to-one correspondence \1/ between the sets LA^lD*) and SSA{D). 
Under the bijection the subset LT^D*) of ZiAy^D*) corresponds with the subset LA{D) of 
SSA{D). Moreover, if {■,-}d* is a linear almost Poisson structure on D* then the corresponding 
skew- symmetric algebroid structure ([■, ■\di Pd) on D is characterized by the following conditions 

\X^\^ = -{X,Y}d*. PD{X){f)orn^ = {forD^,X}D* (2.2) 
forX.Ye V{td) and f G C°^{Q). 

Proof. Let {■, ■}£>* be a linear almost Poisson structure on D* . Then, it is easy to prove that |-, ■\d 



(defined as in (2.2)) is a M-bilinear skew-symmetric bracket. Moreover, since {^-jo* satisfies the 
Leibniz rule, it follows that Pd{X) is a vector field on Q, for X G T{tj:i). In addition, using again 
that {■, -jz)* satisfies the Leibniz rule and Proposition 2.2[ we deduce that 



PD{gX) = gpoiX), for g G C^iQ) and X G T{td). 
Thus, pd '■ ^{td) — ^ X{Q) is a morphism of C°°((5) -modules. 



On the other hand, from (2.2), we obtain that 

{XJYj^ = -{X, if o rz,0^}D* = {pD{X){f) o Tn^)Y - (/ o td^){X, Y}n^. 

Therefore, 

lX,fY]o = flX,Yjn + PDiX)if)Y 
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and (|-, -iDyPn) is an skew-symmetric algebroid structure on td : D ^ Q. It is clear that if {■, -Id* 
is a Poisson bracket then |-, -Jd satisfies the Jacobi identity. 

Conversely, if ([■, -Id, Pd) is an skew-symmetric algebroid structure on to D ^ Q and x E Q 
then one may prove that there exists an open subset U of Q and a unique linear almost Poisson 
structure on the vector bundle r^-i^^^y : t^^{U)* — > U such that 

and 

for X,Y sections of the vector bundle t^^{U) — > U and f,g G C°^{U). Here, (|-, -J^-i^^), p^-i^^-j) is 

the skew-symmetric algebroid structure on t^^{U) induced, in a natural way, by the skew-symmetric 
algebroid structure (|-, ■]d,Pd) on D. In addition, we have that if |-, -Jd satisfies the Jacobi identity 
then {•, is a linear Poisson bracket on t^^{U)*. Thus, we deduce that there exists a unique 



linear almost Poisson structure {■, ■}/)* on D* such that (2.2) holds. □ 



Let {■, be a linear almost Poisson structure on D* and ([■, -Jn, Pd) be the corresponding skew- 
symmetric algebroid structure on td '■ D —>■ Q. If (g*) are local coordinates on an open subset U of 



Q and {X^} is a basis of sections of the vector bundle t^^{U) — > U such that A^* is given by (2.1) 
(on T^\U)) then 

d_ 

C^p and are called the local structure functions of the skew-symmetric algebroid structure 
([■) '}d,Pd) with respect to the local coordinates (g*) and the basis {X^}. 

Next, we will see that there exists a one-to-one correspondence between SSA{D) and the set of 
almost differentials on the vector bundle tj:, : D ^ M. 

Definition 2.6. An almost differential on the vector bundle '■ D —>■ Q is a M.-linear map 

: r(AVB,) ^ TiA'^+'ro*), k e {0, . . . ,n - 1} 

such that 

d^{a A(3) =d^aA(3 + (-l)^a A d^(3, for a G r(A%.) and f3 G r(AVB.)- (2.3) 
// [d^y = then d^ is said to be a differential on the vector bundle T£, : D ^ Q . 

Denote by AD{D) (respectively, D{D)) the set of almost differentials (respectively, differentials) 
on the vector bundle Td : D ^ Q. 

Theorem 2.7. There exists a one-to-one correspondence $ between the sets 2>2>A{D) and ATi{D). 
Under the bisection $ the subset LA{D) of SSA{D) corresponds with the subset D[D) of AD{D). 
Moreover, we have: 

(i) // d^ is an almost differential on the vector bundle td : D ^ Q then the corresponding skew- 
symmetric algebroid structure (|-, ■\di Pd) on D is characterized by the following conditions: 

a{lX,Y\n) = d^{a{Y)){X) - d^{a{X)){Y) - {d^a){X,Y), pD{X){f) = (rf^/)(X), (2.4) 

forX.Ye V{td), a G V{td*) and f e C^{Q). 
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(ii) // ([-j-lDiPr)) is an skew-symmetric algebroid structure on the vector bundle td : D ^ Q 
then the corresponding almost differential is defined by 

k 

(rf^a)(Xo,Xi,...,Xfc) = 5^(-l)Vz?(X,)(a(Xo,...,X„...,Xfc)) 

i=o (2.5) 
+ ^(_l)^+^«(|X„X,]B,Xo,Xi,...,X,,...,X„...,Xfc) 

i<j 

for a G r(AVfl,) and Xq, . . . G T{td). 
Proof. Let be an almost differential on tjj : D ^ Q and suppose that pr, and |-, -Jd are given by 



(2.4). Then, using the fact that 

d^'iff) = fd^'f + f'd^'f, for /, /' G C^iQ), 

we deduce that (|-, -JdiPd) is an skew-symmetric algebroid structure on td : D ^ Q. Moreover, it 
is well-known that if {d^y = then |-, -Jd satisfies the Jacobi identity and the pair ([-, ■]z),Pd) is a 
Lie algebroid structure on td : D ^ Q (see, for instance, [231 US])- 

Conversely, if (|-, ■]d,Pd) is an skew-symmetric algebroid structure on To ■ D ^ Q and d^ is the 



operator defined by (2.5) then, it is clear that 



(rf^«)(X,F) = pz5(X)(a(r)) -pz.(r)(a(X)) -alX,yl^, (rf^/)(X) = pz.(X)(/), 

for / G C°°(M), a G T^td*) and X,Y E T^to). In addition, an straightforward computation proves 
that 

d^{a A(3) = d^aA(3+ {-l)^a A d^(3, for a G T{A''td*) and f3 G r(AV^.)- 

Finally, it is well-known that if |-, -Jd satisfies the Jacobi identity then {d^Y = (see, for instance, 
)■ □ 



Let (|-, -Jd, Pd) be a skew-symmetric algebroid structure on the vector bundle td : D ^ Q and 
d^ be the corresponding almost differential. If (g*) are local coordinates on an open subset U of Q 
and {Xa} is a basis of sections of the vector bundle T^^(f/) — ^ U such that and p^ are the local 
structure functions of the skew-symmetric algebroid structure, then 



d^'x' = pIX", d^X^ = -^C^^X" A X", 



for all i and 7. 



From Theorems 2.5 and 2.7, we conclude the following result 



Theorem 2.8. Let T£, : D ^ Q be a vector bundle over a manifold Q and D* be its dual vector 
bundle. Then, there exists a one-to-one correspondence between the set LA7{D*) of linear almost 
Poisson structures on D* , the set SSA{D) of skew- symmetric algebroid structures on td : D ^ Q 
and the set AD{D) of almost differentials on this vector bundle. 
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3. Skew-symmetric algebroids and the orbit theorem 

Let {D,l-, -Jd, pd) be a skew-symmetric algebroid over Q and be the corresponding almost 
differential. 

We can consider the vector space over M 

H'id'') = {fe C^{Q)/d'>f = 0}. 

Note that if D is a Lie algebroid we have that H^{d^) is the Lie algebroid cohomology group 
associated with D. 

On the other hand, it is clear that if Q is connected and D is a transitive skew-symmetric algebroid, 
that is, 

PniDq) = TgQ, for all g G Q 

then 

H\d^) - M. (3.1) 



Condition (3.1) will play an important role in Section 4.1 



Next, we will see that (3.1) holds if the skew-symmetric algebroid is completely nonholonomic with 
connected base space. 

Let D be the generalized distribution on Q whose characteristic space at the point g G Q is 
It is clear that D is finitely generated. Note that the C°°-module T{D) is finitely generated (see 

HZl)- ^ 

Now, denote by Lie°°(D) the smallest Lie subalgebra of X{Q) containing D. Then Lie°°(D) is 
comprised of finite M-linear combinations of vector fields of the form 

[Xk, [-^fc-i, • • • [X2, Xi] . . . ]] 

with k E N, k 0, and Xi, . . . , Xk G X((5) satisfying 

Xi{q) G bq, for all g G Q 

(see m- 

For each q E Q, we will consider the vector subspace Lie^(D) of TgQ given by 

Lie^(D) = {X{q) G TgQ/X G Lie°°(£))}. 

Then, the assignment 

qeQ^Ue^{b)CTgQ 

defines a generalized foliation on Q. The leaf L of this foliation over the point go ^ Q is ihe orbit 
of D over the point go, that is, 

L = 0...0 0f )(go) G Q/ii eR,Xie X{Q) and Xi{q) G Dg, for all g G Q}. 

Here, 0?' is the flow of the vector field Xi at the time ii (for more details, see |2]). 

Definition 3.1. The skew- symmetric algebroid (D, |-, ■]/), p^i) over Q is said to he completely 
nonholonomic if 

Ue'^ipDiD)) = Lie^(D) = TgQ, for all q G Q. 
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Thus, if Q is a connected manifold, it follows that D is completely nonholonomic if and only if the 
orbit of D over any point Qq ^ Q is Q. 



Remark 3.2. (i) Definition 3.1 may be extended for anchored vector bundles. A vector bundle 
td : D ^ Q over Q is said to be anchored if it admits an anchor map, that is, a vector 
bundle morphism pD : D ^ TQ. In such a case, the vector bundle is said to be completely 
nonholonomic if Ue'^{pD{D)) = Lie^(D) = T^Q, for all q e Q. 
(ii) If D is a regular distribution on a manifold Q then the inclusion map id '■ D —>■ TQ is an anchor 
map for the vector bundle td : D ^ Q. Moreover, the anchored vector bundle td : D — *• Q is 
completely nonholonomic if and only if the distribution D is completely nonholonomic in the 
classical sense of Vershik and Gershkovich [H] . In this sense it is formulated in the literature 
the classical Rashev sky- Chow theorem: If Lie^(-D) = TgQ, for all q G Q, then each orbit 
is equal to the whole manifold Q. 



Now, we deduce the following result 

Proposition 3.3. If the skew-symmetric algehroid {D, [-, ■}d,Pd) over Q is completely nonholonomic 
and Q is connected then H^{d^) ~ M. 

Proof. Suppose that / G C°°{Q) and that d^f = 0. Let go be a point of Q. We must prove that 

v{f) = 0, for all V G Tg^Q. 

The condition {d^f){qo) = implies that v{f) = 0, for all v G Dg^. 
Thus, first we have that 

= MMf)) - MMf)) = [Xi,X2]if), 

for Xi,X2 G X{Q) and Xi{q),X2{q) G Dg, for all qeQ. 

Secondly, since D is completely nonholonomic then Lie^(D) = Tg^^Q. Therefore, there exists a 
finite sequence of vector fields on Q, Xi, . . . such that G Dg, for alH G {1, . . . , k} and 

g e Q, 

= [Xfc,[Xfe„i,...[X2,Xi]...]](go). 
From both considerations, we deduce the result. □ 

However, the condition H^{d^) ~ M does not imply, in general, that the skew-symmetric algebroid 
D is completely nonholonomic. 

In fact, let D be the tangent bundle to 

d d 

If (x,y) are the standard coordinates on M^, it follows that \Xi = —,X2 = 7—} is a global basis 

ox oy 

of r(T]R^) = X(]R^). So, we can consider the skew-symmetric algebroid structure (|-, ■Jtm.'^ , prm^) on 
TM^ which is characterized by the following conditions 

d d 

|Xi,X21tM2 = 0, pTR^iXi) = — , PtK2(^2) = 
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Then, the generahzed distribution D = TM? on is generated by the vector fields 

d ~ d 
ox oy 

Thus, the Lie subalgebra Lie°^(D) of TM^ is generated by the vector fields 

~ d ~ d ~ d 

ox ay ay 

This implies that 

Lie(:o,.o)(^)7^^(^o,.o)K^ ifyo = 0. 

However, if / e C°°(R2) and f = 0, we deduce that 

df df 

— = 0, xy— = 0. 

ox oy 

Consequently, using that / e C°°(R^), we obtain that / is constant. 

Next, we will discuss the case when the generahzed foliation Lie°°{D) ^ TQ. In fact, we will prove 
the following result. 

Theorem 3.4. Let {D. |-, -JdjPd) be a skew- symmetric algebroid over a manifold Q and f be a real 
C^-function on Q such that f = 0. Suppose that L is an orbit of D and that Dl is the vector 
bundle over L given by Dl = Uq^LDq = t^^{L). Then: 

(i) The couple H-, -Jd, Pd) induces a skew- symmetric algebroid structure H', -Jdl, Pd^) on the 
vector bundle td^^ '■ Dl — L and the skew- symmetric algebroid {Dl, {■, -inL, Pd^) is completely 
nonholonomic. 

(ii) The restriction of f to L is constant. 

Proof. It is clear that 

PD{Dq) = C Lie^(L>) = TgL, for all q E L. 
Thus, we have a vector bundle morphism 

PD^ : Dl ^ TL. 

More precisely, pi?^ = {pd)\Dl- 

On the other hand, we may define a R-bilinear skew-symmetric bracket 

l;-]D,:r{DL)xr{DL)^r{DL). 

In fact, if Xl, Yl e ^Dl) then 

lXL,YL]DM = lX,Y]n{q), for aU ? e L 
where X, Y are sections oi td ■ D ^ Q such that 

(Xl) \unL, Y\ur]L — \UnL: 

with U an open subset of Q and q E U. 
Note that the condition 

Lie^(L>) = TgL, for all q e L, 

implies that 

ilX,YjD)\vnL^O 
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for X, y G r(Z}), with V an open subset of Q and 

X{q) = 0, for all g e 1/ n L. 

Therefore, the map |-, -Jd^ is well-defined. We remark that if q & V H L and {Xa} is a local basis 
of T{td) in an open susbset W of Q, with q e FT, such that X = f^Xa in ly then, using that 
(/")|ynVKnL = and that p/3(F)(g) G TgL, we deduce that 

lX,r]^(g) = r{q){X^,Y]u{q)-pD{Y){q){nXM = 0. 
Moreover, the couple {{■, -Jol^ Pdl) is a skew-symmetric algebroid structure on the vector bundle 

In addition, it is clear that 

(Dl), = Pd,{{Dl\) = D,, Lie^(52) = Lie^iD) 

for all q E L. Thus, we have that the skew-symmetric algebroid {Dl,1-, -Jdi^, PDl) is completely 
nonholonomic. 

This proves (i). 

On the other hand, it follows that the condition d^f = implies that 

and, since H^{d^^) ^ M (as a consequence of the first part of the theorem), we conclude that / is 
constant on L. □ 

It will be also interesting to characterize under what conditions there exist functions / G C°°{Q) 



such that {d^Y f = 0. Using Equations (2.4) we easily deduce that: 

((rf^)V) iX,Y) = i[pniX),pniY)]-pnlX,Y]n)f 

for all X,Y E r(r£)). Now, consider the generalized distribution D on Q whose characteristic space 
Dg at the point g G Q is: 

Dg = {{[pniX),pniY)]-pnlX,Yjo)iq) / X,Y e r(r^)}. 

It is also clear that D is finitely generated. Denote by Lie°°{D) the smallest Lie subalgebra of X(Q) 
containing D. Observe that Lie°°{D) C Lie°°(Z)). We deduce that {d^Yf = if and only if / is 
constant on any orbit L oi D. Of course if D is completely nonholonomic then the unique functions 
/ G C°°(Q) satisfying {d^Yf = are / = constant, but it has not to be always the case and it may 
useful to find this particular type of functions in concrete examples. 

4. Linear almost Poisson structures and Hamilton- Jacobi equation 

4.1. Linear almost Poisson structures and Hamiltonian systems. In this section, we will 
consider Hamiltonian systems associated with a linear almost Poisson structure on the dual bundle 
D* to a vector bundle and with a Hamiltonian function on D*. Thus, the ingredients of our theory 
are: 

(i) A vector bundle td : D ^ Q of rank n over a manifold Q of dimension m; 

(ii) A linear almost Poisson structure {■, - Id* on D* and 

(iii) A Hamiltonian function h : D* R on D* . 
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The triplet (D, {-, ■}d*, h) is said to be a Hamiltonian system. 

We will denote by A^y* the linear almost Poisson 2- vector associated with {■, - Id*- Then, we may 
introduce the vector field Oi^^* on D* given by 

jc;;^^* = -t{dh)AD*. 

IK^°* is called the Hamiltonian vector field of h with respect to Ad*. The integral curves of 
IK^^* are the solutions of the Hamilton equations for h. 

Now, suppose that (g*) are local coordinates on an open subset U of Q and that {Xa} is a basis 
of the space of sections of the vector bundle t^^(U) U. Denote by {q\Pa) the corresponding local 
coordinates on D* and by C^^ and the local structure functions (with respect to the coordinates 
(q^) and to the basis {Xa}) of the corresponding skew-symmetric algebroid structure on D. Then, 



using (2.1), it follows that 



Therefore, the Hamilton equations are 

dq^_ :dh_ dpa _ , dh dh 
dt ~ dpa ' dt ~ dq^ ^ "^^^ dp0 ^ ■ 

4.2. Hamiltonian systems and Hamilton- Jacobi equation. Let {D, {-, ■}D*,h) be a Hamilton- 
ian system and a : Q ^ D* he a. section of the vector bundle td* : D* —>■ Q. 

If "K^^* is the Hamiltonian vector field of h with respect to {■,-}d*, we may introduce the vector 
field 'Kf^'^* on Q given by 

<r (?) = (T„(,)r^.)(^J"* («(?))), for q G Q. 



From (4.1), it follows that 



<r(?)epD(/^5), forallgeQ, (4.2) 

where ([■, -JdjPd) is the induced skew-symmetric algebroid structure on the vector bundle td : D ^ 
Q. 

Then, the aim of this section is to prove the following result. 

Theorem 4.1. Let {D,{-,-}o*,h) be a Hamiltonian system and a : Q ^ D* be a section of the 
vector bundle To* : D* ^ Q such that d^a = 0. Under these hypotheses, the following conditions are 
equivalent: 

(i) If c : I ^ Q is an integral curve of the vector field "K^^* , that is, 

cit) = (T,(,(i))rB.)(^^°*(«(c(t)))), for all t E /, 

then a o c : I ^ D* is a solution of the Hamilton equations for h. 

(ii) a satisfies the Hamilton- Jacobi equation 

d^{hoa) = 0. 
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Remark 4.2. Let D be the generalized distribution on Q given hj D = 
the smallest Lie subalgebra of X{Q) containing D. Then, using Theorem 
Hamilton- Jacobi equation holds for the section a if and only if the function hoa : Q ^ M.is constant 
on the leaves of the foliation Lie°°(D). o 



3.4 



(D) and Ue°°{D) be 
we deduce that the 



In order to prove Theorem 4.1[ we will need some previous results: 



Proposition 4.3. Let {■, ■}£,* be a linear almost Poisson structure on the dual bundle D* to a vector 
bundle td : D ^ Q , a : Q ^ D* be a section of Td* : D* ^ Q and #a^, : T*D* TD* be the 
vector bundle morphism between T*D* and TD* induced by the linear almost Poisson 2-vector Ad*. 
Then, a is a 1-cocycle with respect to (i.e., d^a = 0) if and only if for every point q of Q the 
subspace ofTa{q)D* 

^a,D(g) = {T,a){pD{D,)) (4.3) 
is Lagrangian with respect to Ad*, that is, 

#Ao. (('C'a,D(g))°) = ^a,D{q), for all q eQ. 

Remark 4.4. li D = TQ and {■, ■}t*q is the canonical Poisson (symplectic) structure on T*Q then 
a is a 1-form on Q, = d^^ is the standard exterior differential on Q and pd = Ptq '■ TQ — >■ TQ is 



the identity map. Thus, if we apply Proposition 4.3 we obtain that a is a closed 1-form if and only 



if a{Q) is a Lagrangian submanifold of T*Q. This is a well-known result in the literature (see, for 
instance, [Tj). o 

Proposition 4.5. Under the same hypotheses as in Proposition \4-^ if the section a is a 1-cocycle 
with respect to d^ , i.e., d^a = then we have that 

irer(#A,, (a(g))) C (£,,z5(g))°, for all q G Q. 
The proofs of Propositions |4.3| and |4.5| may be found in the Appendix of this paper. 



Proof of Theorem 4-T It is clear that condition (i) in Theorem 4.1 is equivalent to the following fact: 



(i') The vector fields iK^^* and on Q and D* , respectively, are a-related, that is, 

{T,a){%^^: (g)) = :K^* (a(g)), for all q E Q. (4.4) 
Therefore, we must prove that 

(V) ^ (n) 



(ii) Let g be a point of Q. Then, using (4.2), (4.3) and (4.4), we deduce that 

?{^-*(«(g))G£„,z5(g). 



Consequently, from Proposition 4^, we obtain that 

^h"* = #Ad* iVa{q)), for some r]^(^g) E {La,D{q))°- 

Thus, since J-Cf^^* {a{q)) = —ij^{^^,{dh{a{q))), it follows that 

Va(q) +dh{a{q)) E Ker{4fA^,{a{q))). 



Now, using Proposition 4.5 and the fact that rja^q) E (^^^/^(g))^, we conclude that 

dh{a{q)) E (i:a,D(g))°. (4.5) 
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Finally, if aq E Dg, we have that 

d^{h o a){q){ag) = dh{a{q)){{Tqa){pD{aq))) 



which implies that (see (4.3) and (4.5)) 

{h o a){q){aq) = 0. 

(ii) =^ (i') Let g be a point of Q. Then, using that d^{h o a){q) = 0, we deduce that 

dhiaiq)) e (i:a,D(g))°. 

Therefore, it follows that 

^t^'Hq)) = -i^AnMh{a{q))) G #a,. ((£„,D(g))°) 



and, from Proposition 4.3, we obtain that there exists Vq G poi^Dq) C TqQ such that 

%t-*{a{q)) = {Tqa){vq). 

This implies that 

= (r„(,)rz,0(:K^*(«(g)) = Vq 

and, thus, 

□ 

Let (D, {•, ■}D*ih) be a Hamiltonian system and a : Q ^ -D* be a section of the vector bundle 
td. ■.D*-^Q. 

Suppose that (g*) are local coordinates on an open subset U of Q and that {X^} is a basis of 
sections of the vector bundle Tjj^(U) U. Denote by {q\P'y) the corresponding local coordinates on 
D* and by pl^, C^^^ the local structure functions of the skew-symmetric algebroid structure (|-, ■]/), Pd) 
with respect to the local coordinates (g*) and to the basis {X^}. If the local expression of a is 

a{q') = {q\a^{q')) 

then 
and 



A/.oa)=0^p;(g)(|^ +^ 1^ ) = 0, V7, VgGt/. 



Corollary 4.6. Under the same hypotheses as in Theorem 4-i if, additionally, H^{d'^) c:^ or if 
the skew- symmetric algebroid {D, [^JdjPd) is completely nonholonomic and Q is connected, then the 
following conditions are equivalent: 

(i) If c : I Q is an integral curve of the vector field 'K^'^* , that is, 

m = (T„(,(i))rBO(^^°*(«(c(t)))), for all t e /, 

then a o c : I ^ D* is a solution of the Hamilton equations for h. 

(ii) a satisfies the following relation 

ho a = constant . 
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Note that if Q is connected then 
ho a = constant <^= 



dh 



+ 



da^ dh 



0, Vz and Vg G U. 



4.3. Linear almost Poisson morphisms and Hamilton-Jacobi equation. Suppose that td : 
D ^ Q and Tjj : D ^ Q are vector bundles over Q and Q, respectively, and that {■, ■}d* (respectively, 
{tIdO is a linear almost Poisson structure on D* (respectively, D*). Denote by (|-,-]z),Pd) and 
(respectively, (|-, ■]5,P5) and d^) the corresponding skew-symmetric algebroid structure and the 
almost differential on the vector bundle td : D ^ Q (respectively, Tjj : D ^ Q). 

Definition 4.7. A vector bundle morphism (F, F) between the vector bundles td* : D* ^ Q and 



D* 



D* 



Q 



Q 



is said to be a linear almost Poisson morphism if 

for^,ipe C°°{D*). 



(4.6) 



Let {F, F) be a vector bundle morphism between the vector bundles t^* : D* Q and r^. : 
D* Q. If X is a section of Tjj : D —>■ Q then we may define the section {F, F)*X of td '■ D ^ Q 
characterized by the following condition 

aMF, PyXM) = F(a,)(X(F(g))), (4.7) 

for all g G Q and aq & D*. 

Theorem 4.8. Let [F, F) be a vector bundle morphism between the vector bundles td* : D* —>■ Q 
and T[)* : D* ^ Q. Then, {F, F) is a linear almost Poisson morphism if and only if 

l{F,FyX,iF,FyYjn = (F, F)*|X, Fj^, (4.8) 
(TF o p^)((F, FyX) = ps(X) o F, (4.9) 

forX,Y eVirs). 

Proof. Su ppos e that {F, F) is a linear almost Poisson morphism and that Z is a section of : D ^ 
Q. From (4.7), it follows that 



{F,FyZ = ZoF. 



(4.10) 



Now, if X,Y e T{td) then, using (2.2) and (4.10), we deduce that 



|(F, pyx, (F, FyYjD = -{X o F, F o F}^. 
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Thus, from (2.2) and (4.6), we obtain that 



which imphes that (4.8) holds 



On the other hand, if / G C°^{Q) then, using again (2.2) and (4.6), it follows that 



Therefore, from (2.2) and (4.10), we have that 



and, consequently. 



(PS(^)(/) O F) O TD* = pD{{F,FyX){foF) O TD* 

PDiX)if ) o F = pdHF, FyX){foF). 



This implies that (4.9) holds. 



Conversely, assume that (4.8) and (4.9) hold 



Then, if f,g ^ C°°{Q) it is clear that the real functions 

f oTi), o F = f o F otd', g ot^, o F = g o F otd* 
are basic functions with respect to the projection r^* : D* ^ Q. Therefore, we deduce that 
= {(/o r^,) oF,{go T^*) o F}d* = {for^^^go td,}^, o F. 



(4.11) 



Now, if X G T{t£)) then, using (2.2) and (4.10), we obtain that 

{(/ oT5.)oF,Xo F}n* = {{TF o p^){{F, o 



Consequently, from (2.2) and (4.9), it follows that 



{(/ o rsO o F,X o F}d* = {/ o r5,,X}5. o F. 



On the other hand, if "K G r(rj5) then, using (2.2), (4.8) and (4.10), we deduce that 



{XoF,YoF}n* = {X,Y}5,oF. 



Thus, (4.11), (4.12) and (4.13) imply that {F,F) is a linear almost Poisson morphism. 



(4.12) 

(4.13) 
□ 



Let (F, F) be a vector bundle morphism between the vector bundles td* : D* ^ Q and r^* : 
D* — >■ Q. Denote by K^F : A^D* A^Z)* the vector bundle morphism (over F) between the vector 
bundles A'^td* : A'^D* Q and A'^tq* : A'^D* Q induced by F. Then, a section a G T{A^td*) is 
said to be (F, F)-related with a section d G T{A''t£,,) if 

A^F oa = ao F. 

Now, assume that F is a surjective map and that (F, F) is a fiberwise injective vector bundle 
morphism, that is. 



Fq = Fid* ■■ Dl F)^(g) 



is a monomorphism of vector spaces, for all g G Q, and 

F{q) = F{q')^F,{Dl) = F,iD;,). 

Then, we may consider the vector subbundle F{D*) (over Q) of r/), : D* ^ Q. Moreover, if a is a 
section of this vector subbundle we have that there exists a unique section a of To* : D* ^ Q such 
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that a is (F, F)-related with a. In fact, if \oLi\ is a local basis of sections of t^^,-) : F{D*) Q, it 
follows that {a-j} is a local basis of r(r£)*). 

Theorem 4.9. Let {F, F) be a vector bundle morphism between the vector bundles r^* : D* ^ Q 
and TQ* : D* —> Q. 

(i) // {F, F) is a linear almost Poisson morphism then the following condition (C ) holds: 

(C) For each a G T^A^td*) which is {F , F) -related with a G T{A''tq*) we have that 
d^a G r(A'^+VD.) is also {F, F)-related with d^a G r(A^+V5.)- 

(ii) Conversely, if condition (C ) holds, F is a surjective map and {F, F) is a fiberwise infective 
vector bundle morphism then {F, F) is a linear almost Poisson morphism. 

Proof, (i) Suppose that X and Y are sections of : D ^ Q. 



Then, if / G using ( |4.9[ ), we deduce that 

{d^f){X) o F = {poiiF, o F). 



Thus, from (4.7), it follows that 

id^mX) o F =< o o F > . 

Therefore, we have that 

d^foF = Fod^ifoF). 
Now, let a be a section of td* : D* Q which is {F, F)-related with a G r(r£),), that is, 

F o a = a o F. 



(4.14) 
(4.15) 



Then, using (4.8), (4.9) and (4.15), we obtain that 

{d^a){X,Y)oF = pD{{F,FrX){a{{F,FrY))-pD{{F,FrY){a{{F,FrX)) 



which implies that 
This proves that 



-a|(F,F)*X,(F,F)*F]z5 
{d^a){X, Y)oF =< K^F o d^a, {X o F,Y o F) > . 
d^aoF = A^F o d^a. 



(4.16) 



Consequently, from (2.3), (4.14) and (4.16), we deduce the result 



(ii) If X G T{Tf)) and / G C°°{Q) then, using condition (C), we have that 

(p5WoF)(/) =rf^(/oF)((F,F)*X) = {po{{F,FrX)){foF). 
This proves that ( |4.9 ) holds. 

Next, suppose that Y G r(r£)) and that a is a section of td* : D* Q which is (F, F)-related 
with a G r(r/)*). 



Then, from (4.9), it follows that 
a|(F, pyx, (F, FTYh = -{d''a){{F, PyX, (F, F)*F) + p5{X){a{Y)) o F - p5{Y){a{X)) o F 
Thus, using condition (C), we deduce that 

a[(F, pyx, (F, pyYjn = -{d^a){X, Y) o F + p5(X)(«(F)) o F - ps(?)(« W) ° F. 
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This implies that (4.8) holds. □ 



Remark 4.10. Let {F,F) be a linear almost Poisson morphism between the vector bundles td* '■ 
D* ^ Q and tq* : D* ^ Q. Moreover, suppose that F is surjective and that (F, F) is a fiberwise 
injective vector bundle morphism. 



(i) From Theorem 4.9, we deduce that the condition H^{d^) ~ M implies that H^{d^) ~ M. In 
general, the converse does not hold. However, if H^{d^) M and / G C°°{Q) is a F-basic 
function such that d^f = then / is constant. 

ii) If F is a surjective submersion with connected fibers, VgF C Z)g = p^^Dg), for all q E Q, and 
d^f = then / is a F-basic function. Here, VF is the vertical bundle to F. 



Now, we will introduce the following definition. 

Definition 4.11. Let {D, {■,■}£,*, h) (respectively, {D, {■, ■}f)*,h)) be a Hamiltonian system and 
(F, F) he a linear almost Poisson morphism between the vector bundles td '■ D* —>■ Q and Tfj : 
D* Q. Then, (F, F) is said to be Hamiltonian if 

hop = h. 

It is clear that if (F, F) is a Hamiltonian morphism then the Hamiltonian vector fields of h and 
h, and , are F-related, that is, 

(T^F)(<-'(/5)) = J{f *(F(/5)), for all (3 G D*. 

This implies that if /x : / ^ is a solution of the Hamilton equations for h then F o p : I ^ D* is 
a solution of the Hamilton equations for h. 



In addition, from Theorem 4.9, we deduce the following result 



Theorem 4.12. Let {D, {■,■}£)*, h) (respectively, {D, h)) be a Hamiltonian system and {F , F) 

be a Hamiltonian morphism between the vector bundles td* : D* Q and r^, : D* ^ Q. Assume 
that the map F is surjective and that (F, F) is a fiberwise injective vector bundle morphism. 

(i) If a is a section of the vector bundle td* : D* ^ Q such that d^a = 0, it satisfies the 
Hamilton- Jacobi equation for h ( respectively, the strongest condition ho a = constant ) and 
it is (F, F) -related with a e r(r^,) then d^a = and a satisfies the Hamilton- Jacobi equation 
for h (respectively, the strongest condition h o a = constant ). 

(ii) If a is a section of the vector subbundle F{D*) of r/)* : D* ^ Q such that d^a = and 
a satisfies the Hamilton- Jacobi equation for h (respectively, the strongest condition h o a = 

constant ) then d^a = and a satisfies the Hamilton- Jacobi equation for h (respectively, the 
strongest condition ho a = constant ), where a is the section of td* : D* ^ Q characterized 
by the condition F o a = a o F . 

5. Applications to nonholonomic Mechanics 

5.1. Unconstrained mechanical systems on a Lie algebroid. Let ta '■ A Q he a Lie algebroid 
over a manifold Q and denote by (|-, ■]A,Pyi) the Lie algebroid structure of A. 
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If S : A Xq A W is a bundle metric on A then the Levi- Civita connection 

: r(r^) X r(r^) ^ r(r^) 

is determined by the formula 

2g(V|y,Z) = pAiX)i9iY,Z))+pAiY)iSiX,Z))-pAiZ)iSiX,Y)) 

+9{x, lz,YU) + 9{Y, lz,xjA) - Siz, lY,XjA) 

for X,Y, Z & r(A). Using the covariant derivative induced by V^, one may introduce the notion of 
a geodesic of as follows. A curve a : I ^ A is admissible if 

^{ta oa) = pAoa. 

An admissible curve a : I Ais said to be a geodesic if V^(j-)(T(t) = 0, for all t G I. 

The geodesies are the integral curves of a vector field i^g on A, the geodesic flow of A, which is 
locally given by 

Here, (g*) are local coordinates on an open subset V of Q, {Xb\ is an orthonormal basis of sections 
of the vector bundle r^^(?7) — >• ?7, (g*, f ^) are the corresponding local coordinates on A and p^, 
are the local structure functions of A. Note that the coefficients ^ of the connection are 

(for more details, see [Hin])- 

The Lagrangian function L : A ^ M of an (unconstrained) mechanical system on A is defined 

by 

L{a) = -9{a,a) - V{TA{a)) = - V{TA{a)), for a e A, 

V : Q being a real C°°-function on Q. In other words, L is the kinetic energy induced by S 

minus the potential energy induced by V. 

Note that if A is the Liouville vector field of A then the Lagrangian energy El = A(L) — L is 
the real C°°-function on A given by 

^9{a,a) + V{TA{a)) = ^\ 

On the other hand, we may consider the section gradc^V of ta '■ A ^ Q characterized by the following 
condition 

SigradsV^X) = (rfV)(X) = Pa{X){V), VX G T{ta). 

Then, the solutions of the Euler- Lagrange equations for L are the integral curves of the vector 
field on A defined by 

^L = ^S- {gradsVy, 

where (grad^V)^ G X{A) is the standard vertical lift of the section gradt^V. The local expression of 
the Euler-Lagrange equations is 

_ J „B - E _ B C _ j ^ 

q — PbV , V — i^EB^ ^ ^dq^' 



Eiia) = ^9{a,a) + V"(r^(a)) = -||a||^ + l^(r^(a)), for a e A. 
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for all i and E (see [HE]). 

Now, we will denote by bg : A — A* the vector bundle isomorphism induced by S and by ^^g : 
A* A the inverse morphism. If a : Q ^ ^* is a section of the vector bundle ta* : A* —>■ Q we also 
consider the vector field ^L,a on Q defined by 

^L,a(g) = (7'#g(a{g))^A)(^L(#g(a(g)))), for g G Q. 

Corollary 5.1. Let a : Q ^ A* be a 1-cocycle of the Lie algebroid A, that is, d"^a = 0. Then, the 
following conditions are equivalent: 

(i) If c : I ^ Q is an integral curve of the vector field S,L,a on Q we have that j^^o ao c : I ^ A 
is a solution of the Euler-Lagrange equations for L. 

(ii) a satisfies the Hamilton- Jacobi equation 

d^{EL o #g o a) = 0, 

that is, the function -||#g o a||g + V on Q is constant on the leaves of the Lie algebroid 
foliation associated with A. 

Proof. The Legendre transformation associated with the Lagrangian function L is the vector bundle 
isomorphism bg : — *• A* between A and A* induced by the bundle metric S (for the definition of the 
Legendre transformation associated with a Lagrangian function on a Lie algebroid, see [25j). Thus, 
if we denote by S* the bundle metric on A* then, the Hamiltonian function Hl = E^ o ^g induced 
by the hyperregular Lagrangian function L is given by 

HLil) = ^g*(7,7) + ViTA^i^)), for 7 G A*. 

Therefore, if A a* is the corresponding linear Poisson 2- vector on A* and CK^^* is the Hamiltonian 
vector field of Hl with respect to A^* , we have that the solutions of the Hamilton equations are the 
integral curves of the vector field J-C^^* . In fact, the vector fields and are bg-related, that is. 



Tbgoa = <^ 



Consequently, if cr : / — > A is a solution of the Euler-Lagrange equations for L then \)go a : I —>■ A* 
is a solution of the Hamilton equations for H^ and, conversely, if 7 : / ^ A* is a solution of the 
Hamilton equations for Hl then #g07:/^y4isa solution of the Euler-Lagrange equations for L 
(for more details, see [25]). 

In addition, since ta* o bg = r^, it follows that 



^Aq) = (T,(,)rAO(<f («(?))) = <1* (9), for q E Q, 



Thus, using Theorem 4.1 (or, alternatively, using Theorem 3.16 in we deduce the result. □ 



Next, we will apply Corollary 5.1 to the particular case when A is the standard Lie algebroid TQ 
and a is a 1-coboundary, that is, a = dS with S : Q ^ M. a real C°°-function on Q. Note that, in 
this case, the bundle metric S on TQ is a Riemannian metric g on Q and that o a = j^g o dS is 
just the gradient vector field of S, gradgS, with respect to g. 
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Corollary 5.2. Let S : Q ^ be a real -function on Q. Then, the following conditions are 
equivalent: 

(i) If c : I ^ Q is an integral curve of the vector field C,L,ds on Q we have that gradgSoc : I ^ A 
is a solution of the Euler-Lagrange equations for L. 

(ii) S satisfies the Hamilton- J acobi equation 

d{EL o gradgS) = 0, 

that is, the function - \\gradgS\\l + V on Q is constant. 
2 ^ 



Remark 5.3. Corollary 5.2 is a consequence of a well-known result (see Theorem 5.2.4 in p.j). o 

Now, let L : y4 ^ M (respectively, L : A ^ M) be the Lagrangian function of an unconstrained 
mechanical system on a Lie algebroid ta '■ A ^ Q (respectively, '■ A ^ Q) and {F, F) be a linear 
Poisson morphism between the Poisson manifolds {A*, {•, ■}a*) and {A*, {-, ■}a*) such that: 

(i) F : Q ^ Q is a. surjective map. 

(ii) For each q E Q, the linear map Fg = F\a* '■ A* — > A*p(^q) satisfies the following conditions 

S*{F,{f3),F,{(3')) = S*{f3,f3'), for (3,(3' e A;, 
F{q) = F{q')^F,{Al) = FAAl,), 

where 9* (respectively, S*) is the bundle metric on A* (respectively. A*). Note that the first 
condition implies that Fg is injective and an isometry. 

(iii) If V : Q ^ R. (respectively, y : Q — M) is the potential energy of the mechanical system on 
A (respectively. A) we have that V o F = V . 

Then, we deduce that {F, F) is a Hamiltonian morphism between the Hamiltonian systems {A, {-, ■}a* , 
Hl) and {A, {■, ■}a*, Hi), where Hl (respectively. Hi) is the Hamiltonian function on A* (respec- 
tively, A*) associated with the Lagrangian function L (respectively, L). 

Moreover, using Theorem 4.12 we conclude that 

Corollary 5.4. (i) If a : Q ^ A* is a 1-cocycle for the Lie algebroid A (d'^a = 0), it satisfies 
the Hamilton- J acobi equation 

d^{EL o #g o a) = (5.1) 

and it is {F , F) -related with a G V{ta*) then d^a = and a is a solution of the Hamilton- 
Jacobi equation 

rf^(E^ o #g o a) = 0. (5.2) 

(ii) If a : Q F[A*) A* is a 1-c ocycle for the Lie algebroid A (d^a = Oj and it satisfies 
the Hamilton- J acobi equation (5.2) then d^a = and a is a solution of the Hamilton Jacobi 
equation (5.1). Here, a : Q A* is the section of ta* : A* Q characterized by the 
condition F o a = a o F . 

A particular example of the above general construction is the following one. 

Let F: Q^Q = Q/Ghea. principal G-bundle. Denote hj (p : G x Q —* Q the free action of G 
on Q and by T0 : G x TQ — > TQ the tangent lift of 0. T0 is a free action of G on TQ. Then, we 
may consider the quotient vector bundle ta = ttq/g '■ ^ = TQ/G ^ Q = Q/G. The sections of this 
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vector bundle may be identified with the vector fields on Q which are G-invariant. Thus, using that a 
G-invariant vector field is F-projectable and that the standard Lie bracket of two G-invariant vector 
fields is also a G-invariant vector field, we can define a Lie algebroid structure ([^-Ia^Pa) on the 
quotient vector bundle = ttq/g '■ ^ = TQ/G Q = Q /G. The resultant Lie algebroid is called 
the Atiyah (gauge) algebroid associated with the principal bundle F : Q —y Q = Q/G (see 

On the other hand, denote by T*0 : G x T*Q T*Q the cotangent lift of the action 0. Then, 
the space of orbits of T*0, T*Q/G, may be identified with the dual bundle A* to A. Under this 
identification, the linear Poisson structure on A* is characterized by the following condition: the 
canonical projection F : A* = T*Q T*Q/G A* is a Poisson morphism, when on A* = T*Q we 
consider the linear Poisson structure induced by the standard Lie algebroid ta = ttq '■ A = TQ — > Q, 
that is, the Poisson structure induced by the canonical symplectic structure of T*Q (an explicit 
description of the linear Poisson structure on A* ~ T*Q/G may be found in [32]). 

Thus, [F, F) is a linear Poisson morphism between A* = T*Q and A* ~ T*Q/G and, in addition, 
F is a fiberwise bijective vector bundle morphism. 

Now, suppose that S = fi' is a G-invariant Riemannian metric on Q and that V : Q —>■ is a G- 
invariant function on Q. Then, we may consider the corresponding mechanical Lagrangian function 
L : A = TQ M. Moreover, it is clear that g and V induce a bundle metric S on A = TQ / G and a 
real function V : Q ^ and, therefore, a mechanical Lagrangian function L : A = TQ/G M. 

On the other hand, we have that for each g G Q the map Fg : A* = T*Q ^*F(q) — {T*Q/G)F{q) 



is a linear isometry. Consequently, using Corollary 5.4, we deduce the following result 



Corollary 5.5. There exists a one-to-one correspondence between the 1-cocycles of the Atiyah alge- 
broid Ta = Ttq/g '■ ^ = TQ/G ^ Q = Q/G which are solutions of the Hamilton- J acobi equation for 
the mechanical Lagrangian function L : A = TQ/G M and the G-invariant closed 1-forms a on Q 

such that the function -||#o o + ^ is constant. 

An explicit example: The Elroy's Beanie. This system is probably the most simple example 
of a dynamical system with a non-Abelian Lie group of symmetries. It consists in two planar rigid 
bodies attached at their centers of mass, moving freely in the plane (see [29] )■ So, the configuration 
space is Q = SE{2) x with coordinates q = (x, y, 9, ip), where the three first coordinates describe 
the position and orientation of the center of mass of the first body and the last one the relative 
orientation between both bodies. The Lagrangian L : TQ ^ M is 

L = im(i;2 + f) + h^e' + ^2(6 + ^j)' - Vi^P) 

where m denotes the mass of the system and Ji and I2 are the inertias of the first and the second 
body, respectively; additionally, we also consider a potential function of the form V{il)). The kinetic 
energy is associated with the Riemannian metric S on Q given by 

g = m{dx'^ + dy^) + {h + I2)d9'^ + hdO (g) dip + hd^p ® dO + hdiP' . 

The system is S'i?(2)-invariant for the action 

$g(g) = {zi + X cos a — y sina, Z2 + X sina + y cos a,a + 6, ip) 
where g = {zi,Z2,a). 
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Let {^1,^2,^3} be the standard basis of 5e(2), 

[^1,61=0, [6,6] = -6 and [6,6] =6- 

The quotient space Q = Q/SE{2) = {SE{2) x S^)/ SE{2) ~ is naturally parameterized by 
the coordinate ip. The Atiyah algebroid TQ / SE{2) ^ Q is identified with the vector bundle: 

Ta '■ A = TS^ X se{2) S^. The canonical basis of sections of is: s ^,6,^2)^3 j- Since the 

metric S is also 5'£'(2)-invariant we obtain a bundle metric S and a S-orthonormal basis of sections: 



V hh yoip h + h J ^va y/m Vh + h 

In the coordinates {■ijj,v^,v'^,v^,v'^) induced by the orthonormal basis of sections, the reduced La- 
grangian is 

L = \ {{v^r + + {vr + (vr) - . 

Additionally, we deduce that 



Ix,,x4a = o, ■ ix2,x4a = o, 

[X2,X4]^ = , ^3, lX3,X4l^ = —===X2. 

V^i + h v-'i + -'2 



Therefore, the non- vanishing structure functions are 



_ / -'2 ^2 _ / -'2 ^3 _ 1 /-'2 _ 1 



Moreover, 



Pa(^i) = \/^^^^, Pa(^2) = 0, Pa(^3)=0, P^(X4) = 0. 



The local expression of the Euler-Lagrange equations for the reduced Lagrangian system L : A 
is: 



Ii±j2dV 



h{h + h) \fh + T2 



h{h+h) Vh^2 



= 0. 
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From the two first equations we obtain the equation: 

h + h dV 



A section a : A*, a{ip) = {ip,ai{'ilj),a2{ip),a3{ip),a4{ilj)), is a 1-cocycle, i.e. d'^a = 0, 

if and only if a2{ip) = 0, a^i^^) = and -^-^ = 0. Therefore, the Hamilton- Jacobi equation 
d^{Ei o #g o a) = is 



Thus, integrating we obtain 
with ki constant. Therefore, 



dijj dijj 



2y(^) + (ai(V'))' = h 



and all the solutions of the Hamilton- Jacobi equation are of the form 
with ^2 constant. 

5.2. Mechanical systems subjected to linear nonholonomic constraints on a Lie algebroid. 

Let ta : A Q he a Lie algebroid over a manifold Q and denote by (|-, -IajPa) the Lie algebroid 
structure on A. 

A mechanical system subjected to linear nonholonomic constraints on A is a pair 
(L, D), where: 

(i) L : A ^ M is a Lagrangian function of mechanical type, that is, 

L{a) = ^S(a, a) - V'(r^(a)), for a e A, 

and 

(ii) D is the total space of a vector subbundle td : D ^ Q oi A. The vector subbundle D is said 
to be the constraint subbundle. 

This kind of systems were considered in [71 [9l [H] . 

We will denote by in '■ D ^ A the canonical inclusion. We also consider the orthogonal decom- 
position A = D (B D-^ and the associated orthogonal projectors P : A ^ D and Q : A ^ D-^. Then, 
the solutions of the dynamical equations for the nonholonomic (constrained) system {L,D) are just 
the integral curves of the vector field ^(l,d) on D defined by 

^(L,D) = TP O^Loio, 



where is the solution of the free dynamics (see Section 5.1) and TP : TA TD is the tangent 
map to the projector P. 

In fact, suppose that (g*) are local coordinates on an open subset U of Q and that {Xb} = {X^/, Xb} 
is a basis of sections of the vector bundle t^^{U) U such that {X^} (respectively, {X;,}) is an 
orthonormal basis of sections of the vector subbundle t^'^^U) U (respectively, r^liU) U). 
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We will denote by {q^,v^) = {q\v'^,v'') the corresponding local coordinates on A. Then, the local 
equations defining the vector subbundle D are 

= 0, for all b. 

Moreover, if and C§q are the local structure functions of A, we have that the local expression of 
the vector field ^(l,d) is 

^iL,n) = P>^|, - {C^yv^ + Pl^)^- (5.3) 
Thus, the dynamical equations for the constrained system (L, D) are 

q' = pl^v\ = -C'^.yv^ - p\ — , v' = Q. (5.4) 

On the other hand, the constrained connection V : ^{ta) x T{ta) r(r4) associated with the 
system (L, D) is given by 

VxY = PiV%Y) + V|gy, for X,Ye r{rA). 

Therefore, if T^c coefficients of we have that 

1 

2' 



— r^i^ — -^{c^^ + cj^ + c^i,), f",^ — 0. 



Consequently, Eqs. (5.4) are just the Lagrange-D 'Alembert equations for the system {L,D) 
considered in [7] (see also [HI E] ) ■ 

Next, we will introduce a linear almost Poisson structure {■, - Id* on D*. 

Denote by {■, -Ia* the linear Poisson bracket on A* induced by the Lie algebroid structure on A. 
Then, 

{(P,^}d* = Woi*D,^ oi*D}A' o P\ (5.5) 
for if^ip E C°°{D*), where i*^-, : A* ^ D* and P* : D* —>■ A* are the dual maps of the monomorphism 
i£) : D —>■ A and the projector P : D —>■ A, respectively. 

It is easy to prove that {■, ■}d* is a linear almost Poisson bracket on D*. Moreover, if {q^iPB) = 
{q\p-y,Pb) are the dual coordinates of {q\ v^) = {q\ v'^, f*) on A* then it is clear that {q\p-y) are local 
coordinates on D* and, in addition, the local expressions of and P* are 

^*D{(i\p-nPb) = iq\Py), P*{q\Py) = {q\Py,o)- (5.6) 



Thus, from (2.1), (5.5) and ( |5.6 ), we have that 



r ,1 ifdy^dij dip d%l) dip dip 

for (^,^ G C^{D*). 

On the other hand, one may introduce a linear Poisson bracket {■, -Ia on A in such a way that 
the vector bundle map bg : A ^ A* is a Poisson isomorphism, when on A* we consider the linear 
Poisson structure {-, ■}a*- Since the local expression of bg is 

^^{q\v^) = {q\v^) 

we deduce that the local expression of the linear Poisson bracket {■, ■}a is 

f-jx_if9ipdip dip dip E d'P dip 

iip^iP^A-PB^Q^.Q^B g^BQq^) ^BCV Q^BQ^C^ 
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for e 

Using the bracket {■,-}a, one may define a linear almost Poisson bracket {■, ■}nh on D as follows. 
If (f and ijj are real C°°-functions on D then 

We have that 

r~ 7i ifd'fdip d'ip ^d^ dip 

V^U = P.(^^ - - C,s-'Q^,g^s- (5-8) 

Thus, a direct computation proves that {■, -jn/j is just the nonholonomic bracket introduced in 



[7]. Note that, using (5.3) and (5.8), we obtain that ^{l,d) is the Hamiltonian vector field of the 
function {El)\d with respect to the nonholonomic bracket {■, ■}nh, i-e., 

^ = ^(L,D)i'fi) = iEL)\D}nh, 

for ^ e C°^{D) (see also [?]). 

Moreover, if S_d is the restriction of the bundle metric S to Z) and bg^ : D ^ D* is the corresponding 



vector bundle isomorphism then, from (5.7) and (5.8), we deduce that 



{^°^SD^'^°hD}nh = {¥',^}d* obg^, for ip,ip e C°^{D*). 

For this reason, {■■,-}d* will also be called the nonholonomic bracket associated with the 
constrained system {L,D). 

We will denote by ([■, -Jd, Pd) (respectively, d^) the corresponding skew-symmetric algebroid struc- 
ture (respectively, almost differential) on the vector bundle T£, : D ^ Q and by ^^g^ : D* — > D the 
inverse morphism of bg^ : D ^ D*. 



Then, from (2.2) and (5.5), it follows that 



lX,YjD = PltDoX,tnoYjA, Pd{X) = pa{idoX), (5.9) 



for X,Y E T(td)- Therefore, using (2.5), we have that 

d^a = A''i*j^{d'^{P* o a)), for a G T{A''td*). (5.10) 

On the other hand, if a : Q ^ -D* is a section of the vector bundle td* : D* ^ Q one may consider 
the vector field ^(L,D),a on Q given by 

^(L,D)a(g) = (7'#s^Kg))rD)(^{L,D)(#Sz,(a(g)))), for qeQ. (5.11) 

Corollary 5.6. Let a : Q ^ D* be a 1-cocycle of the skew-symmetric algebroid {D, |-, ■]d,Pd), that 
is, d^a = 0. Then, the following conditions are equivalent: 

(i) Ifc : I ^ Q is an integral curve of the vector field ^(^l ^0)01 on Q we have that ^g^oaoc : I ^ D 
is a solution of the Lagrange-D'Alembert equations for the constrained system {L,D). 

(ii) a satisfies the nonholonomic Hamilton- Jacobi equation 

rf^((Ei)|DO#g,oa)=0. 

//, additionally, H^{d^) ~ M ('or the skew -symmetric algebroid (D, [-, -Jd, p^) is completely 
nonholonomic and Q is connected) then conditions (i) and (ii) are equivalent to 

(iii) {El)\d ° 4^Sd ° ^ — constant. 
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Proof. Denote by h(^L,D) ■ D* ^ M. the Hamiltonian function on D* given by /i(L,z)) = {El)\d o #g^, 
by Ad* the hnear almost Poisson 2- vector on D* and by 3^^^^*^^ the Hamiltonian vector field of 
h{L,D) with respect to Ao*. Then, the vector fields ^(l,d) and 5^^^^*^^ on D and D*, respectively, are 



bg^-related. Thus, from ( 5.11[ ) and since td* o bg^ = td, it follows that 

^^(X.)a(^) = (^^W^I^OKX^, («(?))) =e(L,D).(g), 

that is, the vector fields ^/^^^'^ja i{L,D)a are equal. 

Moreover, if a : / ^ -D is a curve on Z), we have that a is a solution of the Lagrange-D'Alembert 
equations for the constrained system (L, D) if and only if bg^ o a : / — is a solution of the 
Hamilton equations for h(^L^£,y 

Therefore, using Theorem 4.1, we deduce that conditions (i) and (ii) are equivalent. 

In addition, if H^{d^) ~ M (or if {D, |-, -IdiPd) is completely nonholonomic and Q is connected) 
then, from Corollary 4.6, it follows that conditions (i), (ii) and (iii) are equivalent. □ 

Remark 5.7. Let be the annihilator of D and 3{D^) be the algebraic ideal generated by 
Thus, a section v of the vector bundle A^A* Q belongs to J(i5°) if 

v{q){yi, . . . ,Vk) =Q, for all g G Q and Vi,...,Vk& Dg. 

Now, let Z{to*) be the set defined by 

Z{td*) = {ae T{TD*)/d^a = 0} 

and 2.(r(£)±)o) be the set given by 

Z(r(Bi)o) = {ae T{^D±Y>)/d^a G J(£'°)} 

where (-0"*")° is the annihilator of the orthogonal complement of D and t^x^o : (-0"*")° ^ Q is 
the corresponding vector bundle projection. Then, using (5.10), we deduce that the map 

Z{t£)*) — » 2.(r(£)±)o), a P* o a 
defines a bijection from Z{to*) on Z{T(^£)±^a). In fact, the inverse map is given by 



a ij^ o a. 



On the other hand, if / is a real C°°-function on Q then 

d^f = Q^{d^f){Q)<ZD\ 

o 

Let {L,D) (respectively, {L,D)) be a nonholonomic system on a Lie algebroid : A ^ Q 
(respectively, : A ^ Q) and {F, F) be a linear almost Poisson morphism between the almost 
Poisson manifolds {D*, {■,-}d*) and (/)*, {■,-}d*) such that: 

(i) F : Q — > Q is a surjective map. 

(ii) For each q E Q, the linear map Fg = F\£)* : D* ^*F(q) satisfies the following conditions 

95,iFgi(3),Fg{(3')) = Sd4P,P'). for G D*, 
F{q) = F{q') ^ Fg{Dl) = FgiD*g,), 
where Sd* (respectively, Sd*) is the bundle metric on D* (respectively, D*). 
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(iii) If V : Q ^ (respectively, : Q — ^ M) is the potential energy for the nonholonomic system 
on A (respectively, A) we have that V o F = V . 

Then, we deduce that {F, F) is a Hamiltonian morphism between the Hamiltonian systems {D, {-, ■}[)* , 
h{L,D)) and (-D, {•, -Id*) ^(l,/)))) where h(^L,D) (respectively, hi^ijj-^) is the constrained Hamiltonian 
function on D* (respectively, D*) associated with the nonholonomic system {L,D) (respectively. 



Moreover, using Theorem 4.12 we conclude that 

Corollary 5.8. (i) If a : Q —>■ D* is a 1-cocycle for the skew-symmetric algebroid D (d^a = 0), 
it satisfies the Hamilton- Jacobi equation 

rf^((Ei)|^O#g^O«)=0 (5.12) 

(respectively, the strongest condition {El)\£i o o a = constant ) and it is {F, F)-related 
with a G r{T£)*) then d^a = and a is a solution of the Hamilton- Jacobi equation 

d^{{Ei\r>o4^S^oa) = {) (5.13) 

(respectively, a satisfies the strongest condition {Ei)\j^ o o a = constant ). 

(ii) If a : Q ^ F{D*) C D* is a 1-cocycle for the skew- symmetric algebroid D (d^a = 0) and 
it satisfies the Hamilton- Jacobi equation (5.13) (respectively, the strongest condition {Ei)\[)0 



^g- o a = constant ) then d a = and a is a solution of the Hamilton Jacobi equation 



(5.12) (respectively, a satisfies the strongest condition {El)\£, o jj^c,^ oa= constant ). Here, 



a : Q ^ D* is the section of To* : D* ^ Q characterized by the condition F o a = a o F . 

5.2.1. The particular case A = TQ. Let L : TQ ^ M be a Lagrangian function of mechanical type 
on the standard Lie algebroid ttq : TQ — > Q, that is, 

Liv) = ^g{v,v) - ViTQiv)), for v G TQ, 

where 5^ is a Riemannian metric on Q and \^ : Q — > M is a real C°°-function on Q. Suppose also 
that D is a distribution on Q. Then, the pair {L,D) is a mechanical system subjected to linear 
nonholonomic constraints on the standard Lie algebroid ttq '■ TQ — > Q. 

Note that, in this case, the linear Poisson structure on A* = T*Q is induced by the canonical 
symplectic structure on T*Q. Moreover, the corresponding nonholonomic bracket {■,-}d* on D* 
was considered by several authors or, alternatively, other almost Poisson structures (on D or on 
t'g(-D) <Z A* = T*Q) which are isomorphic to {■,-}d* also were obtained by several authors (see 
[51II91E21I39]). 

Now, denote by #g : T*Q TQ (respectively, #3^-, : D* ^ D) the inverse morphism of the 
musical isomorphism : TQ ^ T*Q (respectively, b^^ : D —>■ D*) induced by the Riemannian 
metric g (respectively, by the restriction of g to D), hj d the standard exterior differential on 
Q (that is, d = d^^ is the differential of the Lie algebroid ttq ■ TQ Q), by ^(l,d) € X{D) the 
solution of the nonholonomic dynamics and by ^{L,D)a ^ ^(Q) its projection on Q, a being a section 



of the vector bundle td* : D* ^ Q (see (5.11)). Using this notation. Corollary 5.6| and Remark 5.7 
we deduce the following result 

Corollary 5.9. Let a : Q ^ D* be a section of the vector bundle Td* : D* ^ Q such that d{P*oa) G 
3(D°). Then, the following conditions are equivalent: 
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(i) Ifc : I ^ Q is an integral curve of the vector field ^(^l^]:))^ on Q we have that i^g^oaoc : I ^ D 
is a solution of the Lagrange-D'Alemhert equations for the constrained system {L,D). 

(ii) d{{EL)\Doi^g,oa)iQ)CD^. 

Remark 5.10. As we know, the Legendre transformation associated with the Lagrangian function 
L : TQ M is the musical isomorphism : TQ T*Q. Moreover, it is clear that X{Q) C D, 
where X is the vector field on Q given by X = 7^^^ o a. Thus, Corollary |5.9| is a consequence of 



some results which were proved in |20j (see Theorem 4.3 in 120]). On the other hand, if H^{d^) ~ M 
(or if Q is connected and the distribution D is completely nonholonomic in the sense of Vershik and 



Gershkovich [H]) then (i) and (ii) in Corollary 5.9 are equivalent to the condition 



(£'l)|d o o a = constant . 
A Hamiltonian version of this last result was proved by Ohsawa and Bloch [31] (see Theorem 3.1 in 

m)- o 

Remark 5.11. Previous approaches. There exists some different attempts in the literature of 
extending the classical Hamilton- Jacobi equation for the case of nonholonomic constraints [TUl E31 ESI 
Hni 113] ) • These attempts were non-effective or very restrictive (and even erroneous), because, 

in many of them, they try to adapt the standard proof of the Hamilton- Jacobi equations for systems 
without constraints, using Hamilton's principle. See [37] for a detailed discussion on the topic. 

To fix ideas, consider a lagrangian system L : TQ — > M. of mechanical type, that is, L{vq) = 
^9{vg,Vg) — V{q), for Vg G TgQ, and nonholonomic constraints determined by a distribution D of Q, 
whose annihilator is = span{/i^ rfg*}. 

The idea of many of these previous approaches consist in looking for a function S : Q — > M called 
the characteristic function which permits characterize the solutions of the nonholonomic problem. 
For it, define first the generalized momenta 

dS , , 
P'=d^i+ ^^^^ ' 

which satisfy the constraint equations Q^^PifJ.'^ = 0. These last conditions univocally determine A;, as 
functions of q and dS/ dq and therefore we find the momenta as functions 

dS 

V^=V^{q\Q-)■ (5.14) 

By inserting these expressions for the generalized momenta in the Hamiltonian of the system, we 
obtain a version of the Hamilton- Jacobi equation (in its time-independent version): 

dS 

H{q\Pi) = H{q\ —) = constant. (5.15) 
However, if we start with a curve c : I ^ Q satisfying the differential equations 

c^(t) = ^(c'(t),^(c(t)) + V5) (5.16) 
in general, it is not true that the curve ■y(t) = {c\t),pi(t)) is a solution of the nonholonomic equations. 

(5.17) 



This is trivially checked since from Equation (|5.15|) we deduce that: 

dH dH 



dq^ ' dpj 



dq'dqJ dq' ^ dq' 
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but, on the other hand, 



Pi 



d 
dt 



'OS 



dq 

dq^dqi dqi^' 'dqJ 



(5.18) 



solution of the nonholonomic equations (that is, pi 



Substituting Equation (5.17) in Equation (5.18), a curve 7(t) = {d{t),pi{t)) satisfying (5.16) is 



iP+Afe/i^) if it verifies the following condition: 



Ah 



~dq 




0, 5q e Dq. 



(5.19) 



It is well-known (see [361 EZ]) that condition (5.19) takes place when the solutions of the nonholo- 
nomic problem are also of variational type. However, nonholonomic dynamics is not, in general, of 
variational kind (see [H EH [26]). Indeed, a relevant difference with the unconstrained mechanical 
systems is that a nonholonomic system is not Hamiltonian in the standard sense since the dynamics 
is obtained from an almost Poisson bracket, that is, a bracket not satisfying the Jacobi identity (see 

ISl[ISl[221En]). o 



An explicit example: The two-wheeled carriage (see [30]). The system has configuration space 
Q = SE{2) X , where SE{2) represents the rigid motions in the plane and the angles of rotation 
of the left and right wheels. We use standard coordinates (x, 6, ipi: "^2) € SE{2) x T^. Imposing the 
constraints of no lateral sliding and no sliding on both wheels, one gets the following nonholonomic 
constraints: 



xsin^ — ycos6 = 0, 
X cos 6' + ?/ sin 6' + + a-i/'i = 0, 
i; cos 6* + ?/ sin 6* — r6' + a'?/^2 = 0, 

where a is the radius of the wheels and r is the half the length of the axle. 

Assuming, for simplicity, that the center of mass of the carriage is situated on the center of the 
axle the Lagrangian is given by: 

L = \mx^ + Imij' + \je' + \ci,l + \ci,l 



where m is the mass of the system, J the moment of inertia when it rotates as a whole about the 
vertical axis passing through the point (x, y) and C the axial moment of inertia. Note that L is the 
kinetic energy associated with the Riemannian metric g on Q given by 



g = m{dx^ + dy"^) + JdO"^ + Cdipl + Cdi)l 
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The constraints induce the distribution D locally spanned by the following (yf-orthonormal vector 
fields 

1 f d d d a 

Xi = — 2r— a— -arcosfc*- arsm9 — 

Ai V (^Wi ay 

I f d d d 

X2 = — {a^{J -TJiir^)— + {a^J + 4:Cr^ + a^mir^)— + ar{2C + mia'^) — 

d d 
-aia^J + 2Cr2) cos 9— - aia^J + 2Cr'^) sin 9— 

ax oy 

where 



Ai = a/ 4Cr^ + a^J + arriir^ 

A2 = V(a2J + 2Cr2)(2C + mia'^){a'^J + ACr^ + d^r'^m^) 

We will denote by (x, 6', f ^, f ^) the local coordinates on D induced by the basis {Xi,X2}. 

In these coordinates, the restriction, L\y) : D 1 — > M, of L to D is: 

The distribution orthogonal to D is generated by 

= {X3 = tan^— - 7r,X4 = ^— + — + 7^77-,^5 = —jr + ^ + 

ox oy riTLi ox o9 Crotpi ami ux oipi oip2 

Moreover, since the standard Lie bracket [Xi,X2] of the vector fields Xi and X2 is orthogonal to 
D, it follows that (see ( [Kq] )) 

|Xi,X21d = 0, pd{Xi)=Xi, pd{X2)=X2, (5.20) 

where (|-, ■}d,Pd) is the skew-symmetric algebroid structure on the vector bundle To : D ^ Q. 

The local expression of the vector field ^(l,_d) is: 

k{L,D) = + r TT— + 



Ai A2 / dipi A2 dip2 

^ ' ar{2C + mia^)v'^ avi\ d / arv^ cos 9 ^ a{a'^ J + 2Cr'^)v'^ cos 9 \ d 



A2 AiJ 89 \ Ai A2 J dx 

arv^sm9 a(a'^ J + 2Cr'^)v'^ sin 9 \ d 
+ 



Ai A2 J dy 

Furthermore, if {X^,^^} is the dual basis of {Xi,X2} and a : Q ^ -D* is a section of the vector 
bundle Td* : D* ^ Q 

a = aiX^ + a2X'^ , with ai, ^2 e C°°(Q) 

then 

a is a 1-cocycle <^==^ Xi{a2) — X2{ai) = 0. 

In particular, taking 

a = KiX^ + 7^2X2, with Ki, /C2 G M 
trivially is satisfied the 1-cocycle condition. 
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In addition, since El = L, we deduce that 

1 
2 

which imphes that 

-2 , T^2 



{El)\£) o o a = + K2 = constant. 



Thus, using Corollary 5.6, we conclude that to integrate the nonholonomic mechanical system (L, D) 
is equivalent to find the integral curves of the vector field on Q = x x M.'^ given by 

i{L,D),a = KiXi + K2X2 . 

which are easily obtained. 

It is also interesting to observe that, in this particular example, 

Lie°°(D) = {a#^ - adip'^ + 2rrf^}° 

and, thus, D is not completely nonholonomic. From Theorem 3.4 it is necessary to restrict the initial 
nonholonomic system to the orbits of D, that in this case are 

Lk = {{x, y, e, ^1, tlj2) e SE{2) X \ a(V'i - + 2re = k, with keR} 

to obtain a completely nonholonomic skew-symmetric algebroid structure on the vector bundle r/)^^ : 
Dl^. — > Lk- Note that on Lk we can use, for instance, coordinates {x,y,ipi,ip2)- 

5.2.2. The particular case A = TQ/G. Let F: Q^Q = Q/Ghea principal G-bundle and 
= ttq/g '■ A = TQ/G Q = Q/G he the Atiyah algebroid associated with the principal bundle 



(see Section 5.1). 

Suppose that g is a G-invariant Riemannian metric on Q, that : Q — M is a G-invariant real 
C°°-function and that D is a G-invariant distribution on Q. Then, we may consider the corresponding 
nonholonomic mechanical system (L, D) on the standard Lie algebroid = t^q '■ A = TQ — > Q. 

Denote by C,(l,d) ^ ^{D) the nonholonomic dynamics for the system {L,D) and by {■,-}d* the 
nonholonomic bracket on D*. 

The Riemannian metric g and the function V : Q ^ induce a bundle metric S on the Atiyah 
algebroid = ttq/g '■ ^ = TQ/G Q = Q/G and a real G°°-function V : Q ^ M. on Q such that 
V o F = V, where F : Q —y Q = Q/G is the canonical projection. Moreover, the space of orbits D 
of the action of G on D is a vector subbundle of the Atiyah algebroid = ttq/g '■ ^ = TQ/G 
Q = Q/G. Thus, we may consider the corresponding nonholonomic mechanical system {L,D) on 
A = TQ/G. 

Let F : A = TQ — > A = TQ/G be the canonical projection. Then, {F, F) is a fiberwise bijective 
morphism of Lie algebroids and F{D) = D. Therefore, using some results in [7] (see Theorem 
4.6 in [7j) we deduce that the vector field ^(l,d) is Fo-projectable on the nonholonomic dynamics 
^(L,D) ^ ^(-D) of the system (L, D). Here, F^ : D D = D/G is the canonical projection. 

On the other hand, if P : A = TQ ^ D and P : A = TQ /G ^ D = D /G are the orthogonal 
projectors then it is clear that 

FdoP = Pop 

which implies that 

FoP* = FdoP*, (5.21) 
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where F : A* = T*Q A* ~ T*Q/G and Fd : D* ^ D* c:^ D*/G are the canonical projections. 
Moreover, if on A* we consider the hnear Poisson structure induced by the Atiyah algebroid = 
ttq/g ■■ A = TQ/G Q = Q/G then, as we know, F : A* = T*Q ^ A* ^ T*Q/G is a Poisson 
morphism. Thus, using this fact, (5.5) and (5.21), we deduce the following result 

Proposition 5.12. The pair [F^^F) is a linear almost Poisson morphism, when on D* and D* we 

consider the almost Poisson structures induced by the nonholonomic brackets {■,-}d* o-nd {■■,■}[)*, 
respectively. 



Note that Proposition 5.12 characterizes the nonholonomic bracket {■, -js*- 

We also note that the linear map (-Fd)^ = {Fd)\d* '■ D* ^*F(q) — i^* /^)F(q) is a linear isometry. 
for all q & Q. Therefore, from Remark 5.7 and Corollary 5.8, it follows 

Corollary 5.13. Let S the set of the 1-cocycles a of the skew-symmetric algebroid Tfj = t^/g '■ D = 
D/G ^ Q = Q/G which are solution of the nonholonomic Hamilton- Jacobi equation 



(i^((Ez)|5o#g^oa) = 

(respectively, which satisfy the strongest condition {Ei)^j-,oj^g^ oa 
a one-to-one correspondence between S and the following sets: 



constant ). Then, there exists 
Q which 



(i) The set of the G-invariant 1-cocycles a of the skew-symmetric algebroid td : D 
are solutions of the nonholonomic Hamilton- Jacobi equation 

(respectively, which satisfy the strongest condition {El)\d ° #3^ o a = constant ). 

(ii) The set of the G-invariant 1-forms 7 : Q — ^ {D'^)^ ^ T*Q on Q which satisfy the following 
conditions 

dj e 3{D°) and d{EL o #3 o 7)(Q) C Z}0 
(respectively, which satisfy the strongest conditions d'-f G 3{D^) and Eioj^gO'j = constant ). 

An explicit example: The snakeboard. 

The snakeboard is a modified version of the traditional skateboard, where the rider uses his own 
momentum, coupled with the constraints, to generate forward motion. The configuration manifold 
is Q = SE{2) X with coordinates (x, ^, V^, 0) (see [lET]). 




The system is described by a Lagrangian 

L(g, q) = \m{x^ + y^) + hj + 2J{)()^ + \u() + i,f + 
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where m is the total mass of the board, J > is the moment of inertia of the board, Jo > is the 
moment of inertia of the rotor of the snakeboard mounted on the body's center of mass and Ji > 
is the moment of inertia of each wheel axles. The distance between the center of the board and the 
wheels is denoted by r. For simplicity, as in pT], we assume that J + Jq + 2Ji = mr^. 

The inertia matrix representing the kinetic energy of the metric g on Q defined by the snakeboard 

is 

g = mdx^ + mdy^ + mr^dO'^ + Jod9 ® dip + J^dip ® d0 + J^dip^ + 2Jid(f . 

Since the wheels are not allowed to slide in the sideways direction, we impose the constraints 

—X sin(6' + (p) + ij cos(6' + (p) — rO cos = 

—X sin(6' — (p) + ij cos(6' — (p) + rO cos = 0. 

To avoid singularities of the distribution defined by the previous constraints we will assume, in the 
sequel, that ^ ±7r/2. 

Define the functions 

a = — r(cos0cos(6' — 0) + cos0cos(6' + 0)) = — 2r cos^ 0cos6' 
b = — r(cos0sin(6' — 0) + cos0sin(6' + 0)) = — 2r cos^ 0sin6' 
c = sin(20). 

The constraint subbundle To '■ D \ — > Q is 

D - III— — a— b— —\ 
^P^ 1 9-0 '90' dx dy ^86) 

The Lagrangian function and the constraint subbundle are left-invariant under the SE{2) action: 
^gil) = {a + X cos 7 — ?/ sin7, /3 + a; sin7 + y cos 7, •y + 9,ip, 0) 

where g = {a,(3,j) e SE{2). 

We have a principal bundle structure F : Q — > Q where Q = {SE{2) x T^)/SE{2) ~ T^, being 

( d d d ] 

its vertical bundle VF = span < 7^ , 7:— , ttt f- We have that 

ox oy oO J 

r d d d 

S = D nVF = span <^ Y3 = a— + b— + c- 



dx dy 06 



and therefore, 



S^nD = span {¥, = -^,¥2 = ^- = span {y, = ^¥2 = ^- _^(tan0)y3 



dtp k j \ 90' dip 2mr'^ 

where k = m{a^ + 6^ + c^r^) = 4mr^(cos^ 0) (away form = ±7r/2) 
Note that if {^1,^2,^3} is the canonical basis of se(2) then 

d Jo sin 



Y 



r (cos 0)^1 + (sine 



dip mr'^ 
yg = -2r(cos2 0)er + (sin20)S 
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where {i = 1, 2, 3) is the left-invariant vector field of SE{2) such that (e) 
identity element ol SE{2). 

Next, we will denote by {Xi, X2, X^} the ^(-orthonormal basis of D given by 

1 d 



^i, e being the 



where /(0) = Jq — 



Xi = 
X2 = 

^3 = 

Jq sin^ (f) 



1 f d Jo sin ( 



-r (cos 0)^1 + (sin 0)^3 



(cos 0)^1 + - (sin 0)^3 



The vector fields {Xi,X2} describe changes in the internal angles 
and while X^ represents the instantaneous rotation when the internal angles are fixed. 
Consider now the corresponding Atiyah algebroid 

TQ/SE{2) ~ (TT^ X T SE{2))/SE{2) — >Q = T^. 

Using the left translations on SE{2), we have that the tangent bundle of SE{2) may be identified 
with the product manifold SE{2) x 5c(2) and therefore the Atiyah algebroid is identified with the 
vector bundle = : ^ = TT^ x 5e(2) — > T^. The canonical basis of : TT^ x sc(2) — > is 
d d 1 

TTT; 777 5 ^3 r ■ The anchor map and the linear bracket of the Lie algebroid : TT^ X5e(2) — )• 

dip ocp J 

is given by 



d 



PA(e^) = 0, i = l,2,3 



d .8,8 
'dip' Ikp' ^^^dp' ^ dp' 
16,^3^ = -6, 16,^31^ = ^1, 

being equal to zero the rest of the fundamental Lie brackets. 

We select the orthonormal basis of sections, X2, X3, X4, X5}, where 

1 d 



X' 



1 



Jo sin I 



y7(0) \8i^ mr^ 



[-r(cos0){i + (sin 0)^3 



X' 



(cos 0)^1 + - (sin 0)^3 



and {X'^,X'^} is an orthonormal basis of sections of the orthogonal complement to D, D^, with 
respect to the induced bundle metric Sa- 

Taking the induced coordinates {ip, 0, v^,v'^, v^, v^, v^) on TT^ x 5e(2) by this basis of sections, we 
deduce that the space of orbits D of the action of SE{2) on D has as local equations, — Q and 
— 0, being a basis of sections of D, {X[, X'2, X'^} . Moreover, in these coordinates the reduced 
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Lagrangian L : TT^ x se(2) — > M is 

Now, we consider the reduced nonholonomic mechanical system {L,D). 
After, some straightforward computations we deduce that 

r^2Jimf{(p) ry/2Jimf{(p) 
Therefore, the non-vanishing structure functions are: 

Oio — Ooi ==, <-^T? 



r^2J,mf{ct>y rV2Jim/(0)' 
Moreover, 

This shows that P£){D) = TgT"^ and then the skew-symmetric algebroid D — > is completely 
nonholonomic. 

The local expression of the vector field C(l,d) is 

d v"^ d Jq cos I -i d Jo cos (f) 12^ 

Let {(X')^, (X')^, {X'Y} be the dual basis of D* . It induces a local coordinate system: (0, ip,Pi,P2,P3) 
on and, therefore, the non-vanishing terms of the nonholonomic bracket are: 



, ^ Jq cos 

{P1,P2}b* = , P3, 

ry/2Jimf{(p) 

, ^ Jo cos 

Now, we study the Hamilton- Jacobi equations for the snakeboard system. A section a : — > D*, 
a = q;i(0, ip)iX'y + a2{(p, ^){X'f + az{(t), is a 1-cocycle {d^a = 0) if and only if: 

1 da2 1 9q;i JqCOS^ 

= , ^ , ^r— H , =a3 (5.22) 

V2JT 90 /7(0) r^2J,mf{(f>) ^ ' 

1 Sag Jo cos 

= , — , =0^2 (5.23) 

V2J1 90 rv/2Jim/(0) ^ ^ 

= (5.24) 
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Finally, since the skew-symmetric algebroid is completely nonholonomic, the Hamilton- Jacobi 
equation is rewritten as 



+ (0:2(0, + (03(0, = constant 



(5.25) 



Now, we will use this equation for studying explicit solutions for the snakeboard, showing the 
availability of our methods for obtaining new insights in nonholonomic dynamics. 

From Equation (5.24) we obtain that 0:3 = 03(0). Then it is clear that also 0:2 = tt2(0)- Assume 
that ai = constant. Therefore, Equations (5.22) and (5.23) are now, in this case, a system of ordinary 
differential equations: 







da'? 



Jo cos ( 



r^mf{(j)) 
da^ Jo cos 
d(f) r^ymf{(j)) 



"3 



02 • 



(5.26) 
(5.27) 



Moreover, observe that all the solutions of these equations automatically satisfy Equation (5.25) 
since 

da2 , da^ 



Ot2- 



+ "3- 



and a\ = constant 



Solving explicitly the system of equations (5.26) and (5.27) we obtain that 

"2(0) 

JqCi 



C^ivTM+^sm. 



r\ m 



03(0) = 

with Ci,C2 arbitrary constants. Therefore, 



sm 



rwrn 



«(0, i,) = (0, ^; v^Co, Ci + 



"^"^^ sin 0, ^^^^ sin 0-6*2 a//(0)) 



r\/m 



is an 1-cocycle of the skew-symmetric algebroid D T^, for all (Co,Ci,C2) G M.^ and moreover 
it satisfies Equations (5.25). Hence, we can use Corollary 5.6 to obtain solutions of the reduced 
snakeboard problem. First, we calculate the integral curves of the vector field 



D)a- 



0(t) 

m 



Civ/7(0(t)) + ^sin0(t) 



J0C2 



r^mf{(l){t)) 



sin0(i(:) 



whose solutions are: 

0(t) = Cot + C3 

m = c^t- 



C2 



Co 



log 



^C2t sin(C3 



V2 cos(Cot + C3) + \ mr^ - Jq sm\Cot + C3 



C4 (if Co = 0) 



C4 (if Co ^ 0) 



Josin2(C3 
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for all constants Cj G M, 1 < i < 5. Now, by a direct application of the nonholonomic equation we 
obtain that 

v\t) = C,y/f{Cot + + ^ sin(Cot + C-s) 



r\/m 



y^^t) = ^ sin(Cot + C3) - C2VfiCot + C3) 
r^/m 



are solutions of the reduced nonholonomic problem. 

6. Conclusions and Future Work 

In this paper we have elucidated the geometrical framework for the Hamilton- Jacobi equation. 
Our formalism is valid for nonholonomic mechanical systems. The basic geometric ingredients are a 
vector bundle, a linear almost Poisson bracket and a Hamiltonian function both on the dual bundle. 
We also have discussed the behavior of the theory under Hamiltonian morphisms and its applicability 
to reduction theory. Some examples are studied in detail and, as a consequence, it is shown the utility 
of our framework to integrate the dynamical equations. However, in this direction more work must 
be done. 

In particular, as a future research, we will study new particular examples, testing candidates for 
solutions of the nonholonomic Hamilton- Jacobi equation of the form a = d^f, for some / G C°°{Q) 
(if there exists) and moreover we will study the complete solutions for the Hamilton- Jacobi equation 
using the groupoid theory. In this line, we will study the construction of numerical integrators via 
Hamilton- Jacobi theory [TH] . We will also discuss the extension of our formalism to time-dependent 
Lagrangian systems subjected to affine constraints in the velocities. It would be interesting to describe 
the Hamilton- Jacobi theory for variational constrained problems, giving a geometric interpretation 
of the Hamilton- Jacobi-Bellman equation for optimal control systems. Finally, extensions to classical 
field theories in the present context could be developed. 



Appendix 

Let {■, ■}£)* be a linear almost Poisson structure on a vector bundle tjj : D ^ Q, (|-, ■]_d,p_d) be 
the corresponding skew-symmetric Lie algebroid structure on D and a : Q ^ D* he a, section of 
td* : D* ^ Q. If g G Q then we may choose local coordinates (q^) = {q\ g") on an open subset U 
of Q, g G U, and a basis of sections {Xa} = {Xi,X^} of the vector bundle t^^{U) U such that 



Suppose that 



p^(X,)(g) = — , pn{X,){q) = 0. (A.l) 
oq \q 



d 

and that the local expression of a in [/ is 



Pz,(X^)=p^^, lXA,XBjD = C^BXc (A.2) 



a(g^) = (g^Mg^)). (A.3) 
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Denote by Kd* the linear almost Poisson 2-vector on D* and by {q^,PA) = {(1^,(1°', Pi, P-y) the corre- 
sponding local coordinates on D*. Then, from (2.1), it follows that 

d 



AD*(a(g)) 



A 



d 



9q' \a{g) dpi 



C^^{q)ac{q)j— A— . 



(A.4) 



Moreover, using (A.l), (A. 2) and (A. 3), we obtain that 

f)fy ■ 

{d^a){q){X,{q),XM) 
{d^a){q){XM),X,{q)) 
(rf^a)(g)(X,(g),X.(g)) 



dq' I, dq^ I, 

da. 



Cf;iq)aA{q), 



-C^M)<^A{q)- 



(A.5) 



On the other hand, let ^a,D{q) be the subspace of Ta(q)D* defined by (4.3). Then, from (A.l) and 
(A.3), we deduce that 

, r d da A d / . X 

(^Q \a(q) Oq \qOpA\a{q) 



which implies that 

('Ca,D(g))° = (W(«(g)),c^P,(«(g)) 



da^ 



dq' 



dq\a{q)),dp^{a{q)) 



Oq \q 



In addition, using (A.4), one may prove that 

#A,,(dg'^(a(g))) = 



dan 



0, 



#A^.(rfp,(a(g)) - ^ dq\a{q))) 

oq \q 



dq^ 



d ,da 



' -CfAqW{q)) 



-cfJq)Mq) 



dQ' \q 

d 



d_ 

9Pi \a{q) 



^Ao*idpj{a{q)) 



,da 



dp. 



l\a{q) 



' -CfMac{q)) 



dt I, 
-C^,iq)ac{q) 



_d_ 
dp^ 



d_ 

9Pi \a{q) 



Proof of Proposition 4-3. From (A.5), (A.6), (A. 7) and (A. 8), we deduce the result. 
Proof of Proposition \4.5 Suppose that 

Paig) = \udq^{a{q)) + n^dpA{a{q)) E T*,^D*. 



Then, using (A.4) and (A.5), it follows that 



---^ u^, for all i. 

I, 



Thus, from (A.7), we conclude that 



(A.7) 



(A. 



□ 



□ 
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